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I. INTRODUCTION 

In noncommutative quantum field theory, the use of 
coherent states has allowed the implementation of com¬ 
mutation relations among coordinates without violating 
Lorentz invariance 00; further developments have es¬ 
tablished a deep relationship with the noncommutative 
Voros product [3|. The main lesson learned from these 
works is that point-like sources should be replaced by 
smeared distributions, actually Gaussian distributions, 
depending on the noncommutative parameter 9; in the 
limit of vanishing 9 , the standard expressions for Green’s 
functions in quantum field theory are recovered. 

These ideas have been applied recently to gravitational 
systems 0H and the analysis of classical black holes, 
such as the Schwarzschild or the Reissner-Nordstrom 
(RN) spacetime, has been performed. As it is well known, 
these two black holes are special cases of more general 
models involving nonlinear electrodynamics that reduce 
to Maxwell’s electrodynamics in certain limits. 

Among them, Born-Infeld (BI) electrodynamics [T0|, 
EH has received a lot of attention since its formulation, 
arising in several contexts such as in the analysis of the 
low-energy effective action for an open superstring [ 121 
or in the low-energy dynamics of D-branes fl3j ; different 
aspects of the associated black holes have been studied 
in the literature ITdl42Ti | . Additionally, in the presence of 
a gravitational field, the Einstein field equations admit 
a black hole solution [25|; in the static case, the EBIon 
is the corresponding spherically symmetric solution fiE . 
These black hole solutions are the nonlinear generaliza¬ 
tion of the RN black hole solutions; they are described 
by the mass M and the charge Q of the black hole in 
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addition to the BI parameter b, directly related to the 
strength of the electromagnetic field at the position of 
the charge. 

In this paper we construct two noncommutative mod¬ 
els of anti-de Sitter-Einstein-Born-Infeld (adSEBI) space- 
time using smeared distributions for the mass and charge 
of the black hole; the sources of the gravitational field 
have then a noncommutative origin. The interest on 
these models lies first on the natural generalization of 
well-known results on black hole thermodynamics using 
Maxwell’s electrodynamics. Additionally, this approach 
allows us to introduce a fundamental length scale that 
can be related to quantum gravity effects and black holes 
should certainly be affected by them Hill- 

The models provide also a nice illustration of the 
changing behavior of thermodynamical quantities when 
the noncommutative parameter is turn on; the influence 
on critical points is clearly seen and phase transitions also 
exist. Moreover, explicit calculations can be performed in 
spite of the nontrivial nature of both the nonlinear elec¬ 
trodynamics and the noncommutative deformation; the 
amount of work involved in obtaining a solution to the 
noncommutative gravitational field equations is similar 
to that of the usual case. This represents an important 
advantage for the study of quantum effects on gravity 
since we have then a nonpertubative treatment of black 
holes in general that allows close comparison with the 
standard results. 

The paper is organised as follows: In Sec. El we re¬ 
view the Einstein-Born-Infcld (EBI) solution and formu¬ 
late two noncommutative models, first by considering a 
noncommutative distribution of mass and then by using 
a noncommutative distribution of charge; the analysis of 
the metric is presented for both cases. The thermody¬ 
namical properties of each of these models are then ana¬ 
lyzed in Sec. IIIII the equation of state and Gibbs energy 
function for the (3+l)-dimensional case are also found. 
We end up with some remarks and future work in the 
Conclusions. Throughout this paper we will use units 
where G = c = 1. 
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II. NONCOMMUTATIVE MODIFICATIONS OF 
THE FIELD EQUATIONS 

The noncommutative field equations for the model 
can be written down straightforwardly. First, we re¬ 
call that the action of a nonlinear electrodynamics cou¬ 
pled to gravity with cosmological constant A in (d + 1)- 
dimensions is given by 

S = J d d+1 xV~g [R-2A+ L(F )], (1) 

where R is the scalar curvature calculated from the met¬ 
ric g^g := det^—g^) and L(F) is an arbitrary function 
of the electromagnetic tensor F llv . Since Maxwell’s elec¬ 
trodynamics should be valid in the weak field limit we 
must have 


L(F) = -F^F^ + 0{F 4 ) 


( 2 ) 


to lowest order. In the specific case of BI electrodynam¬ 
ics, the function L(F) has the form 


L(F) = 4ir 1- Wl + 


F^Fp 

2b 2 


( 3 ) 


where b is the BI parameter. The classical EBI field and 
conservation equations derived from the above action are 
then (22l 


p 


" V y/1 + F 2 /2b 2 ) 
dL(F) 

R~ ^9^ + Q^ v 


= 0, 


where 


1 


dL(F) 




9 F F a 

crv 1 - i± 


dg^ y/i + F 2 /2b 2 


( 4 ) 


( 5 ) 


Let us consider a static spherically symmetric space- 
time in (d + l)-dimensions 


ds 2 = -e 2 ^ r) dt 2 + e 2v{ - r) dr 2 + r 2 ds^_ 1} , (6) 


where 

ds 2 d _ 1 ) = d 6 \ + sin 2 9\d9\ -\ -+ sin 2 61 ■ ■ • sin 2 9 d _ 2 d9 2 l _ 1 

d -1 ( 3 -1 \ 

= n (n sin2 9i ) de F ( ? ) 

3=1 \i=1 J 

is the line element on the 5 ,d_1 sphere with total area (or 
(d — l)-dimensional solid angle) 


n 


2n d / 2 

(d_1) " re¬ 


using the orthonormal basis 


d° = e^ r) dt, 9 1 = e^ r) dr, d a+1 = r sin 6U d9 a 




( 9 ) 

with a = 1,... ,d— 1, the Ricci components are found to 
be 


Roo = e" 2 " (//'+// 2 + /xV + —/Z 


Ru = e - 2 l '(V'-M' 2 +MV + — 


R a+la+1 = ^^{l-e- 2v ) + - r {-g'+ v')e~ 2v .{lQ) 
We now formulate two noncommutative models. 


A. Model A: Noncommutative mass deformation 


At this level, noncommutativity is introduced in the 
model by defining a noncommutative energy-matter ten¬ 
sor 


:= diag{hi, hi, h 3 ,..., h 3 ), h 3 := (r 2 /ii),r/2r, 

(U) 

where 


hi := -p m {r) = - 


M 


(47 x9) d / 2 


0 -r 2 /i9 


( 12 ) 


and a comma denotes partial derivative with respect to 
the radial coordinate r; the normalization of p m (r) is 
such that 

J d d x p m (r) = M. (13) 

It follows that 

< Z flv =diag(-hi,hi,h 3 ,...,h 3 ), (14) 

in the orthonormal basis given previously. The func¬ 
tion p m in the previous expressions represents a non¬ 
commutative diffused distribution of matter which in the 
commutative limit 0—^0 becomes a point-like source of 
mass M located at the origin. It can be verified straight¬ 
forwardly that this form of the energy-momentum tensor 
satisfies = 0 due to the definition of the function 

h 3 in terms of hi. 

The field equations we are interested in solving are 


R u 


A 


dL{F) 
dg^ 




+ 87T 




(15) 


( 8 ) 


for the purely electric case where F tr = E{r) in the co¬ 
ordinate basis. At this stage the conservation laws in 
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Eqs. d4]) are not being changed from their classical ex¬ 
pressions; this will be done in Sec. Ill Bl by means of a 
noncommutative electric current. 

Let us proceed to the solution for the static electric 
charged case. From the BI Lagrangian we have [22j] 


dL(F) dL(F) 
~d^~ ~ dg 11 ’ 


dL{F) 

dg 22 


dL(F) 
dg 33 


= 0. (16) 


Since the conservation law is the same as in the classical 
case, we have 


L(F) = 4b 2 (l - sJl-Eirf/b 2 ) , (19) 


These expressions are a direct consequence of the form where 
of the BI Lagrangian. From them it follows 

Roo = —-Rn =£* g + v = 0. (If) 

In consequence, we only need to determine one metric 
function. From the field equations for the aa-components 
we obtain 


E(r) 


Q 

^J r 2(d-l) + Q2/ b 2’ 


( 20 ) 


d — 2 1 

„2 ~d -1 V 


d— 2^— 2v\ 


8L[F) ^ 00 




DO 9 ~ 8 nhi. 
(18) 


is the electric field in the coordinate basis due to a charge 
Q located at the origin. Therefore, we arrive to the fol¬ 
lowing differential equation for the function v 


V - ‘’"XT'-' = A - ^ # + + 


or equivalently 


87 tM 
(47 lQ) d / 2 


87 tM 


o -r 2 /40 


(r d ~ 2 e~ 2v ), r = (d- 2)r d ~ 3 - Ar^" 1 + 26(6r d " 1 - VQ 2 + 6 2 r 2(d - 1} ) - 


Straightforward integration gives then 

6 0 1 


e ' 2 " = 1 + ^2 - ^ ^ ( 1 - V 1 + J)2 r 2(d-l) 


( 4 nd) d / 

Q 2 A 2(d — 1) Q 2 


r d-l e -r 2 /ie_ 


d r 


d—2 


ds 


8nM 1 {dr 2 

U7t7 


Jr i/s 2 ( d -i) + Q 2 /b 2 27T d / 2 r d “ 2 V 2 ’ 40 y ’ 

where 7 ( 71 , z) is the lower gamma function 


7 (n,z) := / 

Jo 


= I dtt n ~ 1 e~ t , 


and 60 is an arbitrary constant. Using the identity 

'1 


2 F 1 ( -,c;c+l;-z] =c 


1 J.C—1 


t c ~ l dt 


where c is a constant, the equivalent expression 

bo l.o 2 b 2 r 2 


e -2 " = 1 


r.d—2 


— ; Ar +—— I 1 — \/1 T 


0 (1 + ^) 1/2 ’ 

Q 2 \ 2 (d-l) Q 2 


„ ,'l d-2 3d — 4 Q 2 

x 2Fi -, 


b 2 r 2 ( d ~ 1 '> ' d(d- 2 ) r 2 ( d “ 2 ) 


87tM 1 {dr 2 

777 1 77 I i 


2’ 2(d — 1)’ 2(d — 1)’ 5 2 r 2 ( d_1 )/ 2n d / 2 r d ~ 2 1 \2’49 

I- 


( 21 ) 


( 22 ) 


(23) 


(24) 


(25) 


(26) 


for the solution is obtained. Usually bo = —2 M, where M is the mass of the black hole; in the following we set 








































4 


bo = 0 and therefore, a black hole regular at the origin 
exists. 


B. Model B: Noncommutative BI electrodynamics 


We use now the fact that Eq. m is still valid for the 
noncommutative model since we are only introducing a 
charge density without modifying the form of the BI La- 
grangian. In consequence g+v = 0 and thus the previous 
equation becomes 


We now proceed to analyze the noncommutative defor¬ 
mation of the BI solution for the electromagnetic field. 
For this purpose, the conservation law in Eq. 0 is mod¬ 
ified to 


pw 


where 


3" = [pq(O,0] := 


s/YfWJW 

Q 


(47 l9) d / 2 


= 3", 


0 ~r 2 /40 


(27) 


(28) 


is the noncommutative 4-vector current; the point-like 
charge has been replaced by a smooth distribution in¬ 
volving the noncommutative parameter 9. 

With the Ansatz F tr = E(r ), the only nontrivial con¬ 
servation law to be considered is 


e - 2 ^+ v ) r d- 1 E ( r ) 
sjl - Eirjyi? 


Q 


(47 T0) d /' 


-e A1+l/ r d_1 e _r2/4e . 


(29) 


d r 


l E(r) 


R _ r d-l p ~r 2 /A8 


sJ\-E(rflb 2 J (4tt 
I t follows that 

H(r 


E (r) = 


y/r^d-i) + H(r) 2 /b 2 


where 


H(r) := 


Q 


(4:ir9) d / 2 J 0 


d-1 -z 2 /46 


dz z a ~ l e 


Q f d r 


2 TT d / 2 " 7 V 2 ’ 49 


■ (30) 


(31) 


(32) 


When coupled with the gravitational field, and without 
performing a noncommutative deformation of the mass, 
we obtain the following differential equation for the met¬ 
ric function v 


„d-3 \d-l , r,, 2 id-i H 2 (r) 


(r d - 2 e~^) t r = (d-2)r d ~ i - tlr d ~ l +2b z \r d ~ l - 


b 2 


+ r 2(d—1) 


(33) 


The general solution can then be written as 


e -2 " = 1 + - -hr 2 + 2b 2 [°° ds 

r a ~- d .L 




or equivalently 


r 


(34) 


^ = + r ds -H(s)s d PQ (s) + (d-l)H 2 (s) 

r d ~ 2 d d l V b 2 r 2 ( d ~d J d r d ~ 2 J r ^ s 2 (d-i) + H 2 (s)/b 2 


where we have set bo = —2M; with this choice the met¬ 
ric coefficient goo will blow up at the origin. It can be 
seen that the well known expression for EBI spacetime is 
recovered in the limit 0 —» 0. 


III. THERMODYNAMICAL PROPERTIES 

For each of the black holes solutions in the previous 
section, it is straightforward to obtain their thermody¬ 


namical properties. The first step is to calculate the 
energy density surface for r = R = const .; it is given 

by [33] 


<R) :=®-e 0 (R), 


(36) 


where eo(r) is a energy density surface that can be conve¬ 
niently chosen, and k is the trace of the extrinsic curva¬ 
ture k a t of the (d— l)-hypersurface defined by the metric 
' Furthermore, by writing the metric ds 2 for 
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r = const, as 

- N 2 dt 2 + o ab (dx a + V a dt){dx b + V b dt ), (37) 

we associate to it some metric functions o ab and a shift 
vector V a . 

Once the energy density surface e(R) is known, the 
corresponding total quasilocal energy can be found from 

E(R) = j S d -^xy/ae(R), (38) 

where the integration is over the (d — l)-dimensional hy¬ 
persurface defined by r 2 <is 2 d _ip and o := det{o a b)', this is 
the thermodynamic internal energy within the boundary 
R. Using the above expressions, the temperature [30] 


on the lapse function has been used to write the mass 
M as a function of r# and R, i.e. M = M(rH,R). In 
the following, we discuss the thermodynamical proper¬ 
ties of the noncommutative models given in the previous 
section. 


A. Model A 

First, we review very briefly the calculation of the total 
quasilocal energy E{R)-, we have 


e(R) = 


1 2 N(R) 
87t R 


e o {R) — 


1 N(R) 
47t R 


e 0 (R), (43) 



1 k h 
2t r N(R) ’ 


the pressure 


(39) 


and in consequence 

E(R) = -RN{R) - 4nR 2 e 0 {R). 


(44) 


p- ( dE ^ \ 

■ {din^R*-1))' 

and the specific heat 


C R := 



(?P) / (®S). 

V or H J / V dr H ) 


The function e 0 (R) is chosen as 

(40) 

implying that e(R) vanishes identically for anti-de Sitter 
spacetime (M = 0 = Q). The thermodynamic surface 
pressure is then according to Eq. ODD 


can be found; here kh is the surface gravity at the horizon 
r H , S := Cl /4 is the entropy of the black hole 
and in the last expression above, it is assumed that the 
condition 


N 2 {r H ) = - e 2 ^ rH) = - e ~ 2v{rH) = 0, (42) 


V = 


d 


d{AnR 2 ) 


(RN(R)) + 


d{R 2 e 0 {R)) 

OR 2 


(46) 


In this model the lapse function N 2 (r) is given by the 
right hand side of Eq. (1261) . 


I 

From the condition N(rjj) = 0, the mass M is written as 


8 t tM = 


2n d / 2 r 


-2 

H 


(d rk\ 

46 J L 


1 - i Ar » + i bV 


H 


1-4 1 + 


Q 2 


b 2 r% d ~ l \ 


+ 


2(d — 1) Q 2 
d{d - 2) r 2 ( d - 2 ) 


(l d — 2 

1 

CO 

Q 2 y 

[2’2(d-l) 

2{d — 1) ’ 

6 2 4 (d - 1} )\ 


(47) 


On the other hand, from Eq. (l22l) evaluated at the horizon r#, we obtain the surface gravity of the black hole as 


kh = 


2 r 


d -2 
H 


{d - 2)4" 3 - A4” 1 + 2 b 2 r 


H 


1-4 1 + 


Q 2 \ 87 tM 

W r 2(d-1) J “ U nd \d/2 'H 
H / 


4-l e -^/40 


(48) 


Let us specialize the above expressions to the (3+l)-dimensional case; we have that the temperature at the hypersurface 
r = R = const. > xh is (see Appendix lAl) 


T(R) = 


3 ft 2 

2 ’ le 


N(R) / 3 4 \ 47rr « 
' Yl 2’ 40y 


1 - A 4 + 26 2 4 1-4/1 + 


Q 2 

b 2 r 4 H 


8nM 

{And) 3 / 2 


r 2 H e-^/de 


(49) 
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where 


N 2 (R) = 1 - -A R 2 + -b 2 R 2 1 - \ 1 + 


Q 2 

b 2 R 4 


4 Q 2 
3 R 2 '- 


1 15 
2’ 4’ 4’ 


Q 2 \ _ 


8t tM 1 / 3 R 

7 


6 2 1? 4 7 2tt 3 / 2 R ' V 2 ’ 46» 


and 


8ttM = 


2tt 3 / 2 


Th 


1 H 


2’ 46» 


1 - -A r% + 


b 2 r 2 


H 


1-4/1 + 


6 2 7 


if 


4Q 2 




1 1 5 


Q 2 


Zr 2 H ^ 1 V2’ 4’ 4’ 6 2 - 4 


H 


(50) 


(51) 


Notice the presence of the incomplete gamma functions 
in the expression for the temperature Eq. (H^l) . This is a 
nontrivial modification with respect to the standard case 
and is due to the noncommutative deformation: the ratio 
of the incomplete gamma functions becomes unity in the 
commutative limit 9 —> 0. 


1. Equation of state 


Let us now proceed to find the equation of state for our 
black hole solution in (3+l)-dimensions; we follow idlit . 
We know that the entropy of the solution is given by 

S = \ A i A := 4717 #, ( 52 ) 

with thermodynamic volume and associated pressure 
An r 1 

^ T r »- p : = -fSA < 53 ) 

respectively. Then, by taking T := we have 


v 2nv 2 


-S' 1 ' 


1 + 


+ 


M 


16 

b 2 v A ) ' (4tt(9) 3 / 2 ' 


0 -v 2 /we 


(54) 

where the specific volume v is defined as v := 2r#. A 
typical P — v diagram for a fixed mass M is shown in 
Fig. □ Notice the influence of noncommutativity (solid 
line) on the locus of critical points. 

Critical points are defined by the conditions 



Ta 


= 0, 



(55) 


In the commutative case these two conditions translate 
into a cubic algebraic equation for the specific volume v. 
In our model it becomes the trascendental equation 


- px 


= ae 168 v» 


16 Q 2 - 






i 

X 2 


16Q 2 



P 



FIG. 1: P—v diagram of noncommutative BI black hole (mass 
deformation). The thick lines correspond to the noncommu¬ 
tative equation of state with 9 = 1 for T = 0.04517, 0.035 
and 0.026885 (top to bottom); the dashed lines are the corre¬ 
sponding commutative cases in the same order. We have set 
Q = 1, M = 10 and b = 0.45. 


where 


3 b 2 b 2 8 t xMb 2 

P ''~ 2 5 Q 2 ’ q: ~2 8 Q 4 ’ a 2 2O Q 4 tt 3 / 2 0 7 /2 • (57) 

and 

x~ 2 := v 4 + , \x\<b/4Q. (58) 

In the limit 9 —> 0, the right hand side of Eq. flSBll van¬ 
ishes. 

Fig. [2] shows the graphic solution to Eq. (1561) when 
Q = 1,8t tM = 1 ,b = 1/4 for 9 = 0.1 and 9 = 0.05 
respectively. In the commutative case, for the value of 
b given above, the cubic polynomial on the left hand 
side of this equation has no real roots in the interval 
|x| < b/AQ = 0.0625; no critical points exist in this 
regime. The role of noncommutativity is to change this 
behavior, leading to critical points at x = 0.057328 and 
x = 0.0588325 for 9 = 0.1, where the dashed line crosses 
the x-axis (left figure), and to two critical points located 
at x = 0.605505 and x = 0.0614104 for 9 = 0.05 (right 
figure). 


x 


b 2 


- 24 9 
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2. Gibbs energy 


The Gibbs energy can be readily be found using the 
relationship 


G = P(M adm -TS)=P(M-TS), (59) 


where M is the mass of the source written in terms of the 
horizon radius rn (see Eq. (ED), T the temperature of 
the BI black hole at rjj , S its entropy and /3 a constant. 
Knowledge of the Gibbs function G will allow us to infer 
the presence of phase transitions. Explicitly we have in 
(3+1 )-dimensions 


G = 




_ T JL _ e -r a H /40 

4(4t r6») 3 / 2 



-a 4 


2 6 2 - 3 


H 



(60) 


FIG. 2: In the commutative case (solid line) there are no 
critical points for b < bo = l/y/8Q such that |a;| < 6/4Q; 
noncommutativity modifies this behavior and critical points 
arise as a consequence (8 = 0.1 (left figure) and 9 = 0.05 
(right figure), dashed line). 


where M, as a function of r#, is given in Eq. (15T1) . In 
Fig. [3] we show the Gibbs energy for different values of 
the parameters b and 9. It can be noticed that noncom¬ 
mutativity changes the location of the points where first 
order phase transitions happen. 


B. Model B 

As in the previous section, we have the expression for 
the thermodynamic internal energy 


where eo(-R) is given by Eq. (1451) . We now follow a similar 
procedure as before: first, the lapse function N 2 {r ) is 
given by the right hand side of Eq. (f35ll and therefore, 
the mass M is written as 


E{R) = -RN(R ) - 4t T R 2 e 0 (R), (61) 


together with the pressure 

d 


V = 


<9(4t ri? 2 ) 


(RN(R)) + 


d(R 2 e 0 (R)) 

dr 2 


(62) 


J 


M = 


r 


d—2 

H 


2 



2 b 2 r 2 H 
d 



H 2 {r H ) 

k2 2(rf—1) 


2 1 [°° ^ -H(s)s d pQ(s) + (d-l)H 2 (s) 

d r d H ~ 2 Jr H 3 yV^-l) + H 2 {s)/b 2 


(63) 


using as before the condition N(ru) = 0. The surface gravity, as deduced from Eq. (1351) . is given by 


1 


2 r 


d—2 

H 


(d - 2 )r d H ~ 3 - A 4- 1 + 2 b 2 r d „ l 



H 2 (r H ) 


'1 + 


b 2 4 d ~ l \ 


(64) 


k h = 
































FIG. 3: The plots show the Gibbs function for model A for different values of the BI parameter b and the noncommutative 
parameter 6 , as a function of tempera ture f or fixed pressure and Q = 1. In Fig. |3(a)| the known behavior for the Reissner- 
Nordstrom case is seen meanwhile Fig. 3(g) shows the typical behavior for the Schwarzschild solution. 


For the (3+l)-dimensional case, we have that the temperature at the hypersurface r = R = const, is 

1 


T = 


1 


N(R) 47t th 


1 - A r 2 H + 2 b 2 r 2 


H 


1- Wl- 


where 


H2(r H )X 

-H(s)s 3 pq(s) + 2 H 2 (s) 


nhr )=i - - Ia i - ,/^<?) + hi r * ,_ 

R 3 3 l V b 2 R 4 I 3 R J R s 4 + H' 2 {s)/b 2 


(65) 


(66) 
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and 



1.9 2 , 

1- Ar 2 H + -b 

3 H 3 



H^jh) 

b 2 ri 


2 J_ f°° ^ -H(s)s 3 p Q (s) + 2H 2 (s) 

3 r H J rH S ^/s 4 + H 2 (s)/b 2 


1. Equation of state 


and q{v) are given by 


(67) 


We now obtain the equation of state for this_solution in 
the (3+l)-dimensional case; we again follow (29]. As the 
temperature we take again T = kh/ 27t ; with the same 
definitions of pressure P = — A/8-7T and specific volume 
v = 2 rn as before, we have 


T _1_ [ 1 6H 2 (v) \ 

v 2nv 2 477 y V + b 2 v 4 J 


( 68 ) 


This expression is similar to the standard result but for 
the function H(v) that plays the role of the charge Q. 
As seen in Fig. U in this model noncommutativity has a 
more visible effect on the classical P — v diagram. 


P 



FIG. 4: P — v diagram of noncommutative BI black hole 
(charge deformation). The thick lines correspond to the 
noncommutative equation of state with 9 = 0.1 for T = 
0.04517, 0.035 and 0.026885 (top to bottom); the dashed lines 
are the corresponding commutative cases. We have set Q = 10 
and b = 0.45. 


Critical points are defined again by the conditions 


dP 

dv 


= 0 , 


d 2 P 
~dP 2 ) 


\ = 0 . 


(69) 


We define a new variable 


-2 4 , !6 H 2 {v) 

x := v + 


b 2 


(70) 


and after a long but straightforward calculation we ob¬ 
tain 


x 3 + p[v)x + q(v) = 0, (71) 

as the equation determining the critical points of the 
model. We have again a generalization of the cubic alge¬ 
braic relation for the variable x. Here the functions p{y) 


s(v)p(v) := — ^-[v 2 H' 2 (v/2) + v 2 H[v/2)H''[v/2) 
-6vH{v/2)H'(v/2) + 6H 2 (v/2)\, 

s(v)q(v) := 

s{v) := [2vH(v/2)H'{v/2) - 4H 2 (v/2)] 2 , (72) 


and it is understood that the variable v in these expres¬ 
sions is defined implicitly in terms of x through the re¬ 
lation Eq. 01). Here H' denotes derivative of H with 
respect to its argument. 

It can be seen that in the limit 9 —> 0 the functions p(v) 
and q(v) reduce to the standard expressions —36 2 /32Q 2 
and 36 2 /256Q 2 respectively. In Fig. Owe have plotted 
the curve given by Eq. (1711) for 9 = 0.1 and 9 = 0.2. We 
see that in this case noncommutativity does not lead to 
the appearance of critical points in the region |x| < b/AQ. 



FIG. 5: In the commutative case (solid line) there are no 
critical points for b < bo = l/y/8Q such that |ir| < 5/4Q; in 
the noncommutative model B this feature is still present: the 
plot shows the behavior for 6 — 0.1 (dashed line) and 9 = 2 
(dotted line). 
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2. Gibbs energy 


to find the Gibbs energy. We have 


We use here again the expression 

G = P(M adm -TS)=P(M-TS), (73) 


J 



r H + -A r% - -b 


2„3 

r H 


1 , H\r H ) 

b 2 ri 


4 r , -H{s)s* Pq {s) + 2H\s) 

3 J rH 3 ^s 4 + H 2 (s)/b 2 


(74) 


In Fig. [6] we show the typical behavior of this function for different values of the parameters; it is very similar to 
that of the Schwarzschild solution even when noncommutati vity is present. _ 


IV. CONCLUSIONS 

We have analyzed the influence of noncommutativity 
in the thermodynamic properties of the BI black hole. As 
a first step, we obtained the generic expressions for the 
mass and the temperature as a function of the horizon 
radius for arbitrary dimensions when a noncommutative 
modification of the mass term is performed; this change 
involves the replacement of the classical point like behav¬ 
ior by a smeared distribution. In (3+l)-dimensions, the 
P—v diagrams are quite similar to the commutative case 
and the critical points are found by solving a nonhomoge- 
neous cubic equation, contrary to the classical case where 
the nonhomogeneous term is absent. Since this equation 
is more involved that in the commutative case, an ex¬ 
plicit expression for the critical values is hard to obtain, 
however it can be plotted showing the existence of critical 
points where classically they are forbidden. The Gibbs 
function when plotted as a function of the temperature 
for fixed pressure shows the existence of critical points as 
well and the commutative case is also reproduced; phase 
transitions may or not exist depending on the values of 
the noncommutative parameter. 

We also considered a noncommutative deformation of 
the source giving rise to the BI electromagnetic field. As 
in the first model, the relevant expressions for the mass 
and temperature were obtained for arbitrary dimensions 
and later specialized for (3+l)-dimensions. The P — v di¬ 
agrams have some similarities with the classical case but 
there are also important differences; they show a more 
smooth behavior. This is also seen in the plots of the 
function defining the critical points of the model, where 
noncommutativity makes this function to take higher val¬ 
ues rapidly when compared with the classical case. In 
this case the Gibbs function shows a behavior similar to 
the Schwarzschild solution for 9 and b varying, neverthe¬ 
less some differences can be appreciated. 

From the two models analyzed in this work, model A 
has the interesting feature of enhancing the presence of 
critical points. Plots of the specific heat Cr h are given 
in Fig. [7] using dimensionless variables defined by the 
replacements r# —> rsrn,Q —> rsQ,b —> b/rs,k —> 


A/r|,0 —> r|0; here rg is a constant with dimensions 
of length that may be identified with the Schwarzschild 
radius of the source. 

These graphs show the transition from large to small 
values of b for both commutative and noncommutative 
cases. We notice first that for b = 100 the first criti¬ 
cal point in the noncommutative case is shifted to the 
right compared with the standard case. Furthermore, for 
b = 0.1, the first critical point disappears if 9 is small 
(commutative limit) while it remains for values of 9 of 
order one. Fig. [8] shows the situation in the rn — A plane 
for b = 100. As b takes small values, the minimum of the 
solid line (9 small) goes to (0, -oo) and there is only one 
point where a horizontal line crosses this curve; the non¬ 
commutative case (dashed line) on the contrary remains 
essentially unchanged when going from large to small val¬ 
ues of b and have instead two intersection points with the 
same horizontal line. The use of the specific heat to study 
phase transitions in gravitational models has been done 
previously in the literature [29j, ItO - ItI . 

The thermodynamics of a RN black hole has been stud¬ 
ied previously [34| and it has been shown that this black 
hole behaves as a van the Wals fluid, however the effects 
of quantum fluctuations of spacetime are not considered. 
A straightforward calculation shows that in our model 

A for example, this result can be recovered in the limit 
2 

b —> oo with jg « 1. The equation of state gets then a 
correcting term M/(An9) 3 ^ 2 , leading to the noncommu¬ 
tative expression 

3 67t<3 2 

Pc = 8 + (4t r6») 3 / 2 ' 

for the universal ratio of the critical exponents of the van 
der Wals fluid. 

The techniques used in this work may in principle be 
applied to other kind of metrics describing for example 
rotating black holes with electric and magnetic charge 
given by a nonlinear electrodynamics [35| . It will be of 
interest to pursue this topic in a future work. 
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(a) 







(d) 


(e) 


(f) 





FIG. 6: The plots show the Gibbs fucntion for model B for different values of the BI parameter b and the noncommutative 
parameter 8, as a function of temperature for fixed pressure and Q = 1. Fig. 6(a) reproduces the known behavior for the 
Reissner-Nordstrom case meanwhile in Fig. 6(g) the typical behavior for the Schwarzschild solution is shown. In this case, for 
b fixed, noncommutativity can eliminate some critical points. 
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Appendix A: Calculation of the temperature 

Here we derive a expression for the temperature of the 
black hole in (3 + l)-dimensions used in the main text; 
its generalisation to higher dimensions is straightforward. 
The metric is assumed to be 

ds 2 = -N 2 (r)dt 2 + f~ 2 (r) + r 2 {d0 2 + sin 2 Ode/) 2 ), (Al) 
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with N(r) = f(r). Then the temperature at r = R is [30j | 



A*R 2 1 dp(R) = R 1 dN 2 (R) 
8nR f(R) d(nr 2 H ) AttN(R) drn 


Let us assume that 

N 2 {r) = l-MW(r) + V{r,Q 1 h), (A3) 

where W and V are arbitrary functions of their argu¬ 
ments. Then the condition N 2 {th) = 0 implies 

M= HTS [1 + nr "’°- A)] - (A4) 


FIG. 7: Specific heat C r[1 as a function of the horizon radius 
for b = 100 (left figure) and b = 0.1 (right figure). The solid 
line (black) is the curve for 9 = 10 -2 and the dashed line 
(red) the curve for 9 = 1. The values A = —0.01, Q = 1 are 
fixed in all cases. 



FIG. 8: Curves on the rn — A plane where the specific heat 
C rH diverges for 9 = 10 -2 (black solid line) and 9=1 (red 
dashed line); the dotted horizontal line corresponds to the 
value A = —0.01 and the other parameters have the values 
6 =100, Q=l. 


In consequence 


, , 1 RW(R) d 1 

T ( R ) = T7777T , 75— 777^ I 1 + V ( r H, Q, A)] ■ 


N(R) A-kth drn W{th) 

On the other hand, we have 


(A5) 


dN 2 {r) 

dr 


= -MW r + V r 


W(r H ) 


[1 + V{rH,Q,h)]W >r + V| r (A6) 


and evaluation at r = rn gives 


dN 2 {r) 


dr 


= W{r H ) 


d 


dr H W(r H ) 


[l + V(r H ,Q,R)]. 


It follows that 


(A7) 


T(R) 


1 RW{R) 1 dN 2 {r) 
N(R) r hW( r h) At: dr 

1 RW(R) kh 
N(R) r H W(r H ) 


r H 


(AS) 


This expression was used to obtain Eq. (1491) in the main 
text. 
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